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Using the calculated electron energy band structure of Tl2Ba2Cu06, we compute the dependence of 
the c axis magnetoresistance on the orientation of the magnetic field for different magnitudes of the 
magnetic field. We explain the known experimental results for the in-plane rotation of the magnetic 
field and predict the shape of the magnetoresistance oscillations for the out-of-plane rotations of the 
magnetic field. We show how the latter oscillations can be utilized to reconstruct the shape of the 
Fermi surface and to study the coherence of interplane electron tunneling. 



PACS Numbers: 74.72.Fq, 72.15.Gd 

In a strong magnetic field, the electrical resistivity of 
a layered metal oscillates when the magnetic field is ro- 
tated between the orientations perpendicular and parallel 
to the layers. This effect, called the angular magnetore- 
sistance oscillations (AMRO), was originally discovered 
in the layered organic conductors of the BEDT-TTF fam- 
ily (see review jjj) and remains largely unknown outside 
of the organic conductors research community. AMRO 
should not be confused with the Shubnikov-de Haas mag- 
netoresistance oscillations, which occur when the magni- 
tude, not the orientation, of the magnetic field is changed. 
Using AMRO it is possible to determine not only the 
area, but also the shape of the Fermi surface of a lay- 
ered metal [Q. It would be very interesting to observe 
AMRO experimentally in the high-T c superconductors, 
which are also layered materials. That would provide in- 
formation about the structure of their Fermi surfaces and 
the coherence of interlayer electron motion. 

In this paper, we calculate the dependence of the c- 
axis (interlayer) electrical resistivity of Thil^CuOg on 
the orientation of a magnetic field and obtain detailed 
AMRO curves for comparison with the past and future 
experiments ||. The greater the value of uj c t, where oj c 
is the cyclotron frequency and r is the electron scatter- 
ing time, the more pronounced AMRO are. We believe 
that Tl2Ba2Cu06 is one of the best candidates for the 
experimental observation of AMRO because a relatively 
high value uj c t = 0.9 has been recently achieved in this 
material using a pulsed mag netic field of 60 T Q. More- 
over, AMRO for the magnetic field rotation within the 
most conducting (a, b) plane have been already observed 
experimentally in this material at the field of 13 T and 
lo c t = 0.16 -=-0.31 |§. 

T^I^CuOg has the body-centered tetragonal crystal 
structure with in-plane lattice spacings a = b = 3.6 A 
and the distance d = 11 A between the Cu02 planes. 
The electron band structure of Tl2Ba2CuOg has been 
calculated in Ref . H . Two energy bands cross the Fermi 



level: the Cu-0 hole band centered at the Brillouin zone 
corner X and the Tl-0 electron band centered at V. The 
Tl-0 Fermi surface is a closed spheroid, which does not 
contribute significantly to the conductivity, so it will be 
ignored in our calculations. The Cu-0 band is the generic 
band of the layered high-T c cuprates. Its Fermi surface 
has the shape of a slightly corrugated cylinder, open 
along the c direction. The electron dispersion law is a 
sum of the in-plane and interplane terms: 



e(k x 



k>y> k z ) — £|| \kxi ky) + e±{k x , ky, k z ), 



(1) 



where e is the electron energy, and k x , k y , and k z are 
the electron wave vectors along the a, b, and c axes. 
The in-plane dispersion law en(k x ,k y ) has been calcu- 
lated numerically in Ref. |j) on a mesh of 32 x 32 points 
and interpolated in between. For the interplane disper- 
sion law e±(k x , k y , k z ), we select the tight-binding one, 
in accordance with the body-centered unit cell: 



= — 8<j_ cos(k x a/2) cos(k y b/2) cos(k z d), 



(2) 



where the interplane electron tunneling amplitude t± is 
much smaller than the Fermi energy Ef- 

The orientation of the magnetic field H is characterized 
by the polar angle 9 relative to the c (z) axis perpendicu- 
lar to the layers and the azimuthal angle (j> relative to the 
a (x) axis parallel to the Cu-0 bonds. Within a semi- 
classical picture, the electron wave vector k changes in 
time t according to the Lorentz equation of motion: 

dk/dt = (e/hc) v x H, where v = de/hdk. (3) 

Here e is the electron charge, c is the speed of light, and 
h is the Planck constant. 

The interplane component of the electrical conductiv- 
ity tensor o~ zz , obtained by solving the linearized Boltz- 
mann equation in the r approximation, is given the 
Shockley-Chambers formula fq: 
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Angular oscillations of p Z2 vs 9 for <f) = 0° at 




different values of lo c t. 
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FIG. 2. Angular oscillations of p zz vs 9 for 
different values of ui c r. 
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= 2e* 



d 3 k(0) 



x / dtv z [k(t)}e- t/T . 



(4) 



In Eq. (|J), / is the Fermi distribution function, r is the 
electron scattering time, and the first integral is taken 
over the electron wave vector k(0) that serves as the 
initial condition for determining k(f) from the Lorentz 
equation of motion 

When 8 ^ 0, electrons circle around the Fermi surface 
along the closed orbits obtained by cutting the Fermi 
surface with the planes perpendicular to the magnetic 
field. The different planes are labeled by the component 
kjj of the electron wave vector parallel to the magnetic 
field. In this case, Eq. (^) can be rewritten as 
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47r 3 r c I — exp (— T/r) 

x / dk H dtv z [k{t)} / dt'v z [k(t-t')]e- t,/T 
Jo Jo 



(5) 



where T is the period of electron motion. 

Because t±/Ep is very small, we neglect v z in Eq. (||), 
which then reduces to the in-plane equations of motion 
for k x (t) and k y (t), while k z (t) is determined from the 
geometrical relation 



k z (t) = K z - k^t) t&nS. 
Here K z — kjj/ cos 8, and 

k$(t) — k x (t) cos <p + k y (t) sin <f> 



(6) 



(7) 



is the projection of the in-plane electron wave vector onto 
the in-plane component of the magnetic field M . 

Substituting v z from Eq. (J2j) into Eq. (||), using Eq. 
ml), and taking the integral over K z , we find: 



16e 2 m c t\d 



1 



n 2 K 4 u c cos 1 - exp (~2n /lu c t cos 0) 
x d( dC'cos{d[MC) -fc^(C-C')]tan0} 



afc :c (C) afcy(C) ak x ((-(') aky(( - (') 

x cos cos — - — cos cos — 

2 2 2 2 

x exp(~C /uj c t cos 8). (8) 

In Eq. (|[), i has been replaced by the dimensionless vari- 
able £ = u; c i, where w c is the in-plane cyclotron frequency 
for 8 = 0: 



2neH 



eH 



ch § dki I v m c c 



(9) 



In Eq. (g), the integral is taken along the Fermi surface in 
the (k x , k y ) plane, and m c is, by definition, the cyclotron 
mass. 

We have integrated the Lorentz equations of motion for 
k x (t) and k y (t) and evaluated the integral (||) numerically 
for different values of the parameter uj c t and different ori- 
entations of the magnetic field. The resultant resistivity 
p zz = l/a zz H is shown in Figs, [i], |[ and ||. For a fixed 
azimuthal angle <f>, p zz displays oscillations as a function 
of the polar angle 6. Figs. |f| and |^ demonstrate that, 
while the amplitude of oscillations grows with increasing 
uj c t, the angles 8 n where p zz has the n-th maximum do 
not depend on uj c t. 

The values of 8 n can be related to the shape of the 
Fermi surface via the following analytical argument pj. 
The maxima in p zz become sharp in the limit lo c t I. 
In this limit, as follows from Eq. (0), a zz cx t(TJ^), where 
(■ • •} denotes the averaging over different electron orbits 
(the integration over fctf), and 
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FIG. 3. Angular oscillations of p zz vs 9 for different cf> at 
uj c t — 4.0. The curve for d> = 40° is the same as for 6 — 45°. 
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FIG. 4. The linear dependence of tan6 n vs the maximum 
number n. The lines are offset vertically for different <j>. 



x sin [K z d — fe^(C) d tan 9] 



(10) 



is the average of v z along a given electron orbit. In the 
limit (d/a)tan0 ^> 1, the argument of sine in Eq. ( |ic| ) 
oscillates very fast as a function of £. Thus the integral 
is dominated by the points where the phase is stationary. 
These are the turning points of the electron trajectory, 
where fc^ achieves the maximal fc™ ax and minimal — £;™ ax 
values. Evaluating the integral (|10| ) asymptotically in the 
vicinity of the stationary points, we find 



v z oc sin(K z d) cos(fc™ ax d tan 9 - tt/4). 



(11) 



When the argument of cosine in Eq. ( |ll| ) is equal to 7r(n— 
1/2) (where n = 1,2,.. .), the average electron velocity 
vanishes: v z = 0, thus o zz — > and p zz — ► oo ||. This 
situation corresponds to the resistivity maxima in Figs. 
|l|, |2|, and || and takes place at the angles 



tan 9 n 



ir(n - 1/4) 



The same condition (|T^) follows from the derivation even 
if r varies along the Fermi surface jlO] . 

Eq. (|lj) shows that tan6*„ increases linearly with the 
maximum number n, and the slope of this dependence 
is determined by fc™ ax . In Fig. |[ we plot p zz vs 9 for 
different <$>. For each curve, we determine the angles of 
the resistivity maxima 9 n and plot tan# n vs n in Fig. 
|]. Using Eq. (|l^), we determine fc™ ax from the slopes of 
the lines in Fig. [| and plot /c™ ax vs <f> in Fig. || as the 
black dots. To reconstruct the Fermi surface, it is nec- 
essary to draw a line through each black dot in Fig. ^ 
perpendicular to the vector kj? ax . The envelope of these 
lines gives the Fermi surface ||J . One can check in Fig. || 
that this procedure indeed reproduces the actual Fermi 
surface of our model B, which is shown by the solid line. 



Thus, measuring AMRO in TkiE^CuOe experimentally 
for every angle <f> and applying the described above pro- 
cedure, one could reconstruct the in-plane shape of the 
Fermi surface in this material. 

Now, let us discuss AMRO of p zz vs 4> for the mag- 
netic field rotation within the (x, y) plane at 9 — 90°. In 
this case, the electron orbits are open along the c axis. 
Neglecting v z in Eq. (||), we find that only k z depends 
on time: k z (t) = k z (0) — (eH/cTi) {v x sm(j) — v y cos<fi)t. 
Neglecting e± in Eq. ([!]), we can write the volume of in- 
tegration in Eq. (Q) as d 3 k(0) — dk z (0) de\\ dki/hv, where 
the differential dki is taken along the Fermi surface in the 
(k x ,k y ) plane. Substituting v z from Eq. (g) into Eq. (Q), 
using k z (t), and taking the integrals over en, k z (0), and 
t, we find: 



cr zz (90°^) 



t\rd f dki cos 2 ^f- cos 2 ^ 



v 1 + [w,(0)t] 



(13) 



(12) In Eq. ph, 
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\v x sin 4> — v y cos <t>\ (14) 



is the frequency of the electron motion across the Bril- 
louin zone in the k z direction. Similar equations were 
obtained in Ref. |llf| . 

Using Eq. (|13[), we have numerically calculated 
Pzz (90°, 0) and plotted it in Fig. || for different values 
of uj c t. One can see that p zz is minimal at <f> = 0° and 
maximal at = 45°. This result is in agreement with 
Fig. H and with the experiment || . This behavior can be 
qualitatively understood in the following way. Accord- 
ing to Eqs. (H) and Q, a zz (90°, 4>) is dominated by 
the regions of the Fermi surface where u z (4>) is minimal, 
that is where v is parallel to H |ll[] . When the magnetic 
field points along the x axis, these regions are relatively 
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FIG. 5. Dots: k^ ax vs </>. Solid line: The Fermi surface of 
Tl 2 Ba 2 CuOB. The normals to the radii defined by the dots 
and drawn through the dots envelope the Fermi surface. 
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FIG. 6. Angular oscillations of p zz vs (f> for 
different values of ui c t. 
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flat and large (see Fig. ||). On the other hand, when the 
magnetic field points at 45° to the x axis, these regions 
are curved and relatively small. Thus, the conductivity 
is higher when H is parallel to x . 

The experiment || was actually performed in the 
regime uj c t « 1. In this case, p zz can expanded in the 
powers of the magnetic field H: p zz 



o (0) + Aa (2) 

Pzz \ t-±Pzz 



Ap z 4 z \ where pi°J is the zero-field resistance, and Ap z 2 J 
and Ap z 4 } are the positive terms proportional to H 2 
and H 4 , respectively. Angular dependence appears only 
in the Api 4 } term, which can be written as Api 4 ) — 
Pzz + p^zz cos(4c/>). Expanding Eq. ( |l3| ) in the powers of 
lo c t and calculating the integrals numerically, we find the 
value 0.505 for the dimensionless ratio p z 4 z p{°J /[Api 2 )] 2 , 
which does not depend on the magnetic field and the scat- 
tering time t. This value agrees with the experimentally 
measured one 0.6 ±0.1 |). We also obtain the value 0.16 
for the dimensionless ratio of the angular-dependent and 
angular-independent terms p z 4 J / pi 4 ) . The experiment ||] 
finds that this ratio varies with temperature between 0.15 
and 0.06. The temperature dependence may be due to 
different temperature dependences of r at the different 
parts of the Fermi surface po| , which is not considered 
in our model. 

Our results are obtained from the Boltzmann equation 
in the lowest order in the interplane electron tunneling 
amplitude t±, which appears only as a prefactor in Eqs. 
(||) and (|i~3|). Equivalent results can be obtained using 
the lowest-order perturbation theory in ij_ and the in- 
plane electron Green functions [ JL3| . AMRO exist both 
when the interplane electron motion is coherent t± > H/t 
or weakly incoherent t± < h/r On the other hand, 
if the interplane tunneling is strongly incoherent, so that 
electrons loose phase memory, and the in-plane electron 
momentum is not conserved, AMRO should not exist. 



The experimental observation of AMRO in Tl2Ba2CuOg 
H indicates that the interplane tunneling in this material 
corresponds to the former case. 
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